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ABSTRACT 

We present the results from exploratory Chandra observations of nine high-redshift (z=4.09^1.51) optically 
selected quasars. These quasars, taken from the Palomar Digital Sky Survey (DPOSS), are among the optically 
brightest and most luminous z > 4 quasars known (Mb«— 28.4 to —30.2). All have been detected by Chandra 
in exposure times of w 5-6 ks, tripling the number of highly luminous quasars (Mb < —28.4) with X-ray detec- 
tions at z > 4. These quasars' average broad-band spectral energy distributions are characterized by steeper (more 
negative) a ox values ((a ox }=— 1.81 ±0.03) than those of lower-luminosity, lower-redshift samples of quasars. We 
confirm the presence of a significant correlation between the ultraviolet magnitude and soft X-ray flux previously 
found for z > 4 quasars. The joint w 2-30 keV rest-frame X-ray spectrum of the nine quasars is well parameter- 
ized by a simple power-law model whose photon index, T » 2.0±0.2, is consistent with those of lower-redshift 
quasars. No evidence for significant amounts of intrinsic absorption has been found (Nh ^ 8.8 x 10 21 cm~ 2 at 
90% confidence). In general, our results show that z ~ 4. 1 — 4.5 quasars and local quasars have reasonably similar 
X-ray and broad-band spectra (once luminosity effects are taken into account), suggesting that the accretion mech- 
anisms in these objects are similar. We also present near-simultaneous optical spectra for these quasars obtained 
with the Hobby-Eberly Telescope; this is the first time optical spectra have been published for seven of these 
objects. The objects presented in this paper are among the best z > 4 targets for X-ray spectroscopy with XMM- 
Newton and next-generation large-area X-ray telescopes. These will detect or constrain iron Ka emission lines 
down to rest-frame equivalent widths of w 50 eV and intrinsic column densities down to Nn ~ a few x 10 21 cm -2 
at z« 4. We also present 45 new ROSAT upper limits for z > 4 quasars and a likely (3a) HRI detection of the 
blazarGB 1713+2148 at z = 4.01. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: general — X-rays: galaxies 



1. INTRODUCTION 

At z ~ 0-2, X-ray emission is thought to be a universal prop- 
erty of quasars, and X-rays have also been studied from many 
z w 2-4 quasars. In contrast to radio (e.g., Schmidt et al. 1995; 
Stern et al. 2000; Carilli et al. 2001a), millimeter (e.g., Omont 
et al. 1996, 2001; Carilli et al. 2001b), sub-millimeter (e.g., 
McMahon et al. 1999; Priddey & McMahon 2001; Isaak et 
al. 2002), and optical (e.g., Schneider, Schmidt, & Gunn 1989; 
Kennefick et al. 1995; Fan et al. 2001; Constantin et al. 2002) 
wavelengths, where many studies of z > 4 quasars have been 
conducted, observations in the X-ray regime are limited but im- 
proving rapidly (see Brandt et al. 2002b for a review). Prior 
to 2000, there were only six published z > 4 X-ray detections. 
These were for a heterogeneous mixture of objects not suit- 
able for statistical studies because of their limited number and 
varying selection criteria. In the last two years, the number of 
z > 4 X-ray detected quasars has increased to ks 40 (e.g., Kaspi, 
Brandt, & Schneider 2000; Vignali et al. 2001, hereafter V01; 
Brandt et al. 2002a,b; Bechtold et al. 2002). 3 

If z > 4 quasars are radiating close to their maximum (Ed- 
dington) luminosities, their black holes require the assembly 
of a mass > 10 8 -10 9 M (e.g., Haiman & Cen 2002) in a 
relativistic potential when the Universe was only « 10% of 
its present age. X-ray studies of high-redshift quasars repre- 



sent a powerful, direct way of revealing the physical conditions 
in the immediate vicinities of the first massive black holes to 
form in the Universe. The comparison of their X-ray properties 
with those of local quasars may shed light on the physical pro- 
cesses by which quasar central power sources evolve over cos- 
mic time. If these quasars are indeed characterized by high ac- 
cretion rates, as suggested by some models (e.g., Kauffmann & 
Haehnelt 2000), they could plausibly provide evidence for some 
energetically important phenomena such as accretion-disk in- 
stabilities (e.g., Lightman & Eardley 1974; Shakura & Sunyaev 
1976; Gammie 1998) and "trapping radius" effects (e.g., Begel- 
man 1978; Rees 1978). 

The comparison of high-redshift quasar optical-to-X-ray 
spectral energy distributions (SEDs; quantified by means of 
0C ox , the slope of a nominal power law connecting 2500 A and 
2 keV in the rest frame) with those of lower-redshift samples 
suggests a steepening of the continuum shapes at high redshifts 
(V01; Bechtold et al. 2002). However, some studies indicate 
that the results may depend on the sample selection criteria and 
on the choice of the X-ray data used in the analysis (e.g., Kaspi 
et al. 2000; Brandt et al. 2002b). Further observations of ho- 
mogeneous samples of quasars are clearly required, especially 
at the highest redshifts. Moreover, by determining the typical 
X-ray fluxes and oc ox values of high-redshift quasars one should 
be able to assess more reliably the fraction of ionizing flux pro- 
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vided by quasars at the epoch of re-ionization (e.g., Becker et 
al. 2001; Djorgovski et al. 2001; Fan et al. 2002; Pentericci 
et al. 2002), although X-rays alone are not expected to produce 
a fully ionized intergalactic medium (e.g., Donahue & Shull 
1987; Venkatesan, Giroux, & Shull 2001). 

In this paper we present the X-ray properties of a sam- 
ple of nine optically bright z > 4 quasars observed with the 
Chandra X-ray Observatory during Cycle 3. These quasars 
were targeted because they are among the most optically lu- 
minous objects in the Universe [see Fig. 1 for a comparison 
with the known z > 4 quasars; note that in the cosmology 
adopted in this paper, the bright nearby (z = 0.158) quasar 
3C 273 would have an absolute B-band magnitude of -26.4]. 
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FIG. 1. — Absolute B-band magnitude versus redshift for the known z > 4 
quasars (i.e., with Mb < —23). Quasars with an X-ray detection or tight upper 
limit are marked with larger symbols. The PSS quasars presented in this paper 
are plotted as filled triangles, while open triangles indicate quasars observed 
by Chandra and published in previous works (e.g., V01; Brandt et al. 2002a). 
The two filled squares indicate z > 4 quasars serendipitously discovered by 
Chandra at z = 4.93 (Silverman et al. 2002) and z = 5.18 (Barger et al. 2002). 
ROSA 7+observed quasars with tight constraints in the Kaspi et al. (2000) sam- 
ple (see also V01) and in Appendix A are plotted as filled circles, while a 
BALQSO at z = 5.74 (Brandt et al. 2001b) and a z = 4.45 X-ray selected 
radio-quiet quasar in the Lockman Hole (Schneider et al. 1998) are plotted as 
open circles. Blazars (e.g., Kaspi et al. 2000; one of the five objects is pub- 
lished here for the first time; see §3.1 and Appendix A for details) are shown 
as open stars. The uncertainties on Mb (typically f= 0.4 magnitudes) are dom- 
inated by systematic errors due to the uncertainties in the optical magnitude 
measurements, the assumed optical continuum slope, and quasar variability. 



The choice of these quasars was also motivated by the possi- 
bility of obtaining near-simultaneous optical observations (use- 
ful to constrain their SEDs) with the Hobby-Eberly Telescope 
(HET). Although the number of known z > 4 quasars has dra- 
matically increased in the past few years (e.g., Djorgovski et 
al. 1998; Anderson et al. 2001; Peroux et al. 2001; Storrie- 
Lombardi et al. 2001), it is unlikely that large numbers of new 
z > 4 quasars will be found with optical luminosities larger than 
those of the quasars presented here. For example, the most lu- 
minous z > 4 quasar discovered by the Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000) to date is SDSS 0244-0816 with 

4 See http://www.astro.caltech.edu/~george/z4.qsos for a listing of high-redshift quasars. 

5 High-quality optical spectra of PSS 0133+0400 and PSS 0209+0517 have been published by Peroux et al. (2001). 



an absolute B-band magnitude of w —29.1. The quasars pre- 
sented in this paper (M B « —28.4 to -30.2) are thus likely 
to comprise a significant fraction of the most luminous z > 4 
optically selected objects. Since optical magnitude is corre- 
lated with the soft (0.5-2 keV) X-ray flux (see, e.g., Fig. 4 
of V01), these objects are ideal targets for snapshot observa- 
tions with Chandra. They are probably also among the best 
z > 4 quasars to select for future X-ray spectroscopic observa- 
tions with XMM-Newton; at present there is only a handful of 
z > 4 quasars known to be sufficiently X-ray bright for effective 
XMM-Newton spectroscopy. Our targets have been selected 
from the Palomar Digital Sky Survey (DPOSS; e.g., Djorgovski 
et al. 1998; hereafter these objects will be referred to as PSS), 4 
which covers the entire northern sky (8 > —3°) and can find 
quasars up to z « 4.6 (PSS 1347+4956 at z = 4.51 is the high- 
est redshift quasar discovered by the DPOSS thus far). They 
were identified via their distinctive red colors and confirmed 
with optical spectroscopy. We also present HET optical photo- 
metric and spectroscopic observations for all of the quasars in 
our sample. These were obtained nearly simultaneously with 
the Chandra observations, in order to provide constraints on 
these quasars' SEDs and moderate-quality optical spectra. This 
is the first time optical spectra have been published for seven of 
these objects. 5 

Throughout this paper we adopt Hq=70 km s -1 Mpc~' in 
a A-cosmology with i2M=0.3 and Q.a=0.7 (e.g., Lineweaver 
2001). For comparison with the values reported in V01 (with 
Hq=70 km s _1 Mpc -1 and Q.m=1), the A-cosmology increases 
the luminosities typically by a factor of « 2.3-2.4 in the red- 
shift range 4-5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra observations and data reduction 

The nine z > 4 optically selected quasars have been observed 
by Chandra during Cycle 3. The observation log is reported in 
Table 1 . All of the sources were observed with the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2002) with 
the S3 CCD at the aimpoint. Faint mode was used for the event 
telemetry format, and ASCA grade 0, 2, 3, 4 and 6 events were 
used in the analysis. No background flares occurred during the 
observations, the background being constant to within 10-20% 
in all observations. The average 0.5-2 keV background value 
is « 7.6 x 10~ 8 counts s -1 arcsec~ 2 . 

Source detection was carried out with WAVDETECT (Dobrzy- 
cki et al. 1999; Freeman et al. 2002). For each image, we cal- 
culated wavelet transforms (using a Mexican hat kernel) with 
wavelet scale sizes of 1, 1.4, 2, 2.8, 4, 5.7, 8, 11.3 and 16 pix- 
els. Those peaks whose probability of being false were less 
than the threshold of 10~ 6 were declared real; detections were 
typically achieved for the smaller wavelet scales of 1 .4 pixels 
or less as expected for these distant sources. Due to the excel- 
lent angular resolution of Chandra (with a point spread function 
FWHM< 1") for on-axis positions, we were able to locate the 
X-ray sources very accurately (see Table 1), thus avoiding pos- 
sible mis-identifications. X-ray positions have been determined 
using S EXTRACTOR (Bertin & Arnouts 1996) on the 0.5-2 keV 
images and are consistent with those obtained from WAVDE- 
TECT. Source searching was performed in four energy ranges: 
the ultrasoft band (0.3-0.5 keV), the soft band (0.5-2 keV), 
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Table 1 

Chandra AND HET OBSERVATION LOG 



Object 






Optical 


Optical 


Aopt-x a 


Chandra 


Chandra b 


HET 


Name 


z 


R 


012000 


82000 


(arcsec) 


Obs. Date 


Exp. (ks) 


Obs. Date 


PSS 0121+0347 


4.13 


18.3 


01 21 26.1 


+03 47 07 


1.0 


2002 Feb 07 


5.69 


2001 Nov 19 


PSS 0133+0400 


4.15 


18.0 


01 33 40.3 


+04 01 00 


0.4 


2001 Nov 26 


6.08 


2001 Nov 25 


PSS 0209+0517 


4.14 


17.6 


02 09 44.6 


+05 17 13 


0.1 


2002 Jan 14 


5.78 


2001 Dec 09 


PSS 0926+3055 


4.19 


16.5 


09 26 36.3 


+30 55 05 


0.7 


2002 Mar 11 


5.29 


2001 Dec 09 


PSS 0955+5940 


4.34 


18.3 


09 55 11.3 


+59 40 31 


0.4 


2002 Apr 14 


5.68 


2001 Dec 19 


PSS 0957+3308 


4.20 


18.0 


09 57 44.5 


+33 08 22 


0.6 


2001 Dec 24 


6.04 


2001 Dec 10/14 


PSS 1326+0743 


4.09 


16.9 


13 26 11.9 


+07 43 58 


0.7 


2002 Jan 10 


5.89 


2001 Dec 19 


PSS 1347+4956 


4.51 


17.4 


13 47 43.3 


+49 56 21 


0.3 


2002 May 02 


5.89 


2002 Feb 07 


PSS 1443+5856 


4.26 


17.5 


14 43 40.7 


+58 56 53 


0.3 


2002 Jun 18 


5.79 


2002 Mar 14 



NOTE. — The optical positions of the quasars have been derived using the source detection package SEXTRACTOR (Bertin 
& Arnouts 1996) on the the Palomar Digital Sky Survey (DPOSS) images. An average of the redshifts of the emission lines 
in HET spectra (excluding Lycc; see Table 4) was used to determine the redshifts of the quasars. Only two quasars have 
published redshifts: PSS 0133+0400 (z = 4.154; Peroux et al. 2001) and PSS 0209+0517 (z = 4.174; Peroux et al. 2001). 
For the latter quasar, we obtain a slightly different redshift because we averaged two emission lines (excluding Lya), while 
Peroux et al. (2001) used Lya only. R-band magnitudes have been obtained from HET photometric observations. 
a Distance between the optical and X-ray positions. 
b The Chandra exposure time is corrected for detector dead time. 



the hard band (2-8 keV), and the full band (0.5-8 keV); at the 
average redshift of z ~ 4.2, these energy bands correspond to 
the 1.6-2.6, 2.6-10.4, 10.4-41.6, and 2.6-41.6 keV rest-frame 
bands, respectively. 

All of the quasars have been detected with exposure times 
of 5.3-6.1 ks. The WAVDETECT photometry measurements are 
shown in Table 2; we have checked these results with manual 
aperture photometry (using a 2" radius circular cell) and found 
good agreement between the two techniques. Using a larger 
(4" radius) aperture gives results consistent with those reported 
in Table 2 within the errors. Table 2 also reports the hard- 
ness ratios [defined as (H — S)/(H + S), where S is the soft-band 
(0.5-2 keV) counts and H is the hard-band (2-8 keV) counts], 
the band ratios (i.e., the ratio of the 2-8 keV to 0.5-2 keV 
counts), and the effective X-ray power-law photon indices (F) 
derived from these band ratios assuming Galactic absorption. 
All of the quasars have effective photon indices consistent, 
within the errors, with those of z ~ 0-2 quasars in the rest- 
frame 2-10 keV band (r w 1.7-2.3; e.g., George et al. 2000; 
Mineo et al. 2000; Reeves & Turner 2000). To check the re- 



sults derived from the band ratios, we also performed rough 
spectral analyses for all of the sources using the Cash statistic 
(Cash 1979) with XSPEC (Version 1 1.2.0; Arnaud et al. 1996) 
in the 0.5-8 keV band. Although the number of counts for each 
source is too low to provide tight constraints and the signal-to- 
noise ratio dramatically decreases at energies above 2 keV, we 
found general agreement in the power-law photon indices de- 
rived with this method and those obtained from band ratios. The 
inclusion of the 0.3-0.5 keV counts provides for some sources 
slightly different photon indices. However, it must be noted 
that Chandra ACIS suffered from calibration uncertainties in 
the ultrasoft band (G. Chartas 2002, private communication) 
even before the 2002 quantum efficiency loss discovery; there- 
fore possible differences must be taken with caution. 

The soft-band images of the nine quasars, adaptively 
smoothed at the 2a level using the algorithm of Ebeling, White, 
& Rangarajan (2002), are shown in Fig. 2. The crosses mark the 
optical positions of the quasars. The differences between the 
optical and X-ray positions of the quasars are generally < 1" 
(see Table 1). 



Table 2 

X-ray Counts, Hardness Ratios, Band Ratios, and Effective Photon Indices 







X-ray Counts" 








Object 


[0.3-0.5 keV] 


[0.5-2 keV] [2-8 keV] 


[0.5-8 keV] 


Hardness Ratio b Band Ratio b 


r b 



PSS 0121+0347 2.0+ii 58.81^ 9.8+^ 68.51^ -0.71l£„ 0.17+^ 2.1 

PSS 0133+0400 3.0±f;| 31.9±^ 3.9±f;| 36.8±£j — 0.78±g;|^ 0.12±g;^g 2.4_ 05 

PSS 0209+0517 <4.7 21.0+^ < 3.0 22.9+|| < -0.75 < 0.14 > 2.3 

PSS 0926+3055 7.0±|-| 45.9+^| 13.9+| f 59.8jf? -0.54+g-^ 0.30±g"J| 1.6+0.3 

^-H 9.0 fi <6.4 I0.<r g 

PSS 0957+3308 < 3.0 18.0+^;| < 3.0 17.9^1 < -0.71 < 0.17 > 2.0 



0.(1 



PSS 0955+5940 5.01, t 9.0 + ? o < 6.4 10.91, ? < -0.17 <0.71 > 0.8 



PSS 1326+0743 8.0+^ 44.9+™ 14.8+^| 60.6±?;| — 0.50jlg;}f 0.33±j];}g 1.5+0.3 

PSS 1347+4956 < 6.4 27.0tf 3 2 4.0+ig 30.9+ff -0.74^ 0.14+gJg 2.2+JJj 

PSS 1443+5856 < 3.0 7 -°-2 6 < 6A 10.0+^ < -0.04 < 0.91 > 0.6 

a Errors on the X-ray counts were computed according to Tables 1 and 2 of Gehrels (1986) and correspond to the lo~ level; these 
were calculated assuming Poisson statistics. The upper limits are at the 95% confidence level and were computed according to 
Kraft, Burrows, & Nousek (1991). 

b Errors on the hardness ratios [defined as (H — S)I(H + 5), where S is the soft-band (0.5-2 keV) counts and H is the hard-band 
(2-8 keV) counts], the band ratios (i.e., the ratio of the 2-8 keV to 0.5-2 keV counts), and the effective photon indices are at 
the ks lo" level and have been computed following the "numerical method" described in § 1.7.3 of Lyons (1991); this avoids the 
failure of the standard approximate variance formula when the number of counts is small (see § 2.4.5 of Eadie et al. 1971). 
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FIG. 2. — Chandra soft-band (0.5-2 keV) images of the nine high-redshift PSS quasars presented in this paper. In each panel, the horizontal axis shows the Right 
Ascension, and the vertical axis shows the Declination (both in J2000 coordinates). Each image is 97" x 97". The images have been adaptively smoothed at the 2c 
level using the algorithm of Ebeling et al. (2002). Crosses mark the optical positions of the quasars. The possible extension of PSS 0926+3055 is discussed in §3.5. 



2.2. Hobby-Eberly Telescope observations 

Optical spectroscopic and photometric observations of all 
nine quasars have been carried out with the 9-m Hobby-Eberly 
Telescope (Ramsey et al. 1998). The data were obtained with 
the HET's Marcario Low-Resolution Spectrograph (LRS; Hill 
et al. 1998a,b; Cobos Duenas et al. 1998) within at most w 22 
days (in the rest frame) of the Chandra observations (see Ta- 
ble 1). High-redshift, high-luminosity quasars (2 < z < 3.4) 
have longer rest-frame optical variability time scales than low- 
redshift, low-luminosity quasars (z < 0.4), and their variations 
over « 1.5 years (in the rest frame) amount to « 10% (Kaspi 
2001). Therefore, we do not expect our targets to have varied 
significantly (in the optical) during the time between the op- 



tical and X-ray observations. The spectroscopic observations 
were taken with the 300 line mm -1 grating, an OG515 block- 
ing filter, and a slit width of 2-3" (depending on the seeing 
conditions), providing spectra from 5100-10200 A at a resolv- 
ing power of « 240-300. The exposure time was chosen to 
be 900 s for all but one of quasars, the exception being for 
PSS 0957+3308 (two close observations of 800 and 900 s, re- 
spectively). The spectroscopic data were reduced using stan- 
dard IRAF routines. Bias exposures taken on each night were 
stacked and subtracted, and bias-subtracted frames were then 
fiat fielded in a standard manner using internal lamp fiats ob- 
tained during the same night. The wavelength calibration was 
performed using a Ne arc lamp, while the flux calibration was 
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FIG. 3. — Optical spectra of the nine PSS quasars presented in this paper obtained with the Hobby-Eberly Telescope at a resolving power of ks 240-300. The 
principal emission lines are indicated. Some of them are not reported in Table 4 because of the uncertainties in their measurements. The arrows in the panel for 
PSS 1443+5856 mark the Si IV and C IV broad absorption lines. 



achieved using one standard star for each night. The calibrated 
spectra are shown in Fig. 3. 6 Given that there are flux losses 
due to the slit size and that we did not demand photometric 
conditions for our optical observations, there is considerable 



uncertainty in the flux density normalization. 

To address the long-term optical variability of the present 
sample of quasars, we also obtained 120 s 7?-band exposures 
immediately before the spectroscopic observations. Each im- 



6 The HET optical spectra are available at http://www.astro.psu.edu/users/niel/papers/papers.html. 
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age covered a field w 4' on a side, and using United State Naval 
Observatory (USNO) photometry of DPOSS objects it was pos- 
sible to obtain an approximate photometric calibration. Two 
objects (PSS 0926+3055 and PSS 0955+5940) showed signif- 
icant variations (see §4). 

3. X-RAY ANALYSIS 

The principal X-ray, optical, and radio properties of our tar- 
gets are given in Table 3. A description is as follows: 
Column (1 ). — The name of the source. 
Column (2). — The Galactic column density (from Dickey & 
Lockman 1990) in units of 10 20 cm" 2 . 

Column (3). — The monochromatic rest-frame AB 14 5o(i +z ) 
magnitude (defined in §3b of Schneider et al. 1989) with 
estimated errors ^0.1 mag, corrected for Galactic red- 
dening. AB 145 o(i +z ) magnitudes have been derived from 
HET /?-band magnitudes assuming the empirical relationship 
AZ?i45o(i +z )=/?— 0.684 z + 3.10 (effective in the redshift range 
covered by our sample); see Table 2 of Kaspi et al. (2000). 
AB 14 5 ( 1+Z ) magnitudes have not been derived directly from the 
spectra because of flux-calibration uncertainties. 
Columns (4) and (5). — The 2500 A rest-frame flux density 
and luminosity. These were computed from the AB 1450(1 + z ) 
magnitude assuming an optical power-law slope of a=— 0.79 
(S v °< V 01 ; Fan et al. 2001) to keep consistency with V01; see 
V01 for further details about the choice of the optical slope. 
Changing the power-law slope of the optical continuum from 
a=-0.79 (Fan et al. 2001) to oc=-0.5 (e.g., Vanden Berk et al. 
2001) reduces the 2500 A rest-frame flux densities and lumi- 
nosities by w 15%. 

Column (6). — The absolute B-band magnitude calculated ac- 
cording to the formula 

M B =ABi 450( i +z) -51ogD L + 5+2.51og(l+z)-0.83 (1) 
where Dl is the luminosity distance. We note that changing 
the power-law slope of the optical continuum from a=— 0.79 
to a=— 0.5 gives fainter absolute B-band magnitudes by w 0.35 
magnitudes. 

Columns (7) and (8). — The observed count rate in the 
0.5-2 keV band and the corresponding flux (corrected for 
Galactic absorption). The 0.5-2 keV counts have been con- 
verted into fluxes (using PIMMS Version 3.2d; Mukai 2002) 
assuming a power-law model with T = 2 (as derived for sam- 
ples of lower redshift quasars; e.g., George et al. 2000; Mineo 
et al. 2000; Reeves & Turner 2000). Changes of the photon in- 
dex in the range F =1.7-2.3 lead to only a few percent change 
in the measured X-ray flux. The soft X-ray flux derived from 
the 0.5-2 keV counts has also been compared with that derived 
using the full-band (0.5-8 keV) counts. This comparison gave 
similar results to within typically 10-20%. 
Columns (9) and (10). — The rest-frame 2 keV flux density and 
luminosity. These were computed assuming the same power- 
law photon index as in the count-rate-to-flux conversion. 
Column (11). — The 2-10 keV rest-frame luminosity corrected 
for Galactic absorption. 

Column (12). — The optical-to-X-ray power-law slope, a ox , 
defined as 

„ _ log(/2kev//2500 A) 
log(V 2k ev/V2500 A) 



where fa k eV and / 2 500 A WQ tne ^ ux densities at 2 keV and 
2500 A, respectively. 7 The reported errors have been computed 
at the k lo level following the "numerical method" described 
in § 1.7.3 of Lyons (1991). Both the statistical uncertainties on 
the X-ray count rates and the effects of possible changes in the 
X-ray (r w 1.7-2.3) and optical (aw -0.5 to -0.9; e.g., Van- 
den Berk et al. 2001; Schneider et al. 2001) continuum shapes 
have been taken into account (see §3 of V01 for further de- 
tails). We note that changing the power-law slope of the optical 
continuum from a=— 0.79 to a=— 0.5 induces a small increase 
(w 0.03) in the a ox values. 

Column (13). — The radio-loudness parameter (e.g., Keller- 
mann et al. 1989), defined as R=fs GHZ//4400 A ( rest frame). 
The rest-frame 5 GHz flux density was computed from the 
FIRST (Becker, White, & Helfand 1995) or NVSS (Condon 
et al. 1998) observed-frame 1.4 GHz flux density assuming 
a radio power-law slope of a = —0.8. The upper limits re- 
ported in the table are at the 3a level. The rest-frame f^oo a 
flux density has been computed from the AB 1450 (i +z ) magnitude 
assuming an optical power-law slope of a = —0.79. Radio- 
loud quasars (RLQs) are characterized by R > 100, whereas 
radio-quiet quasars (RQQs) have R < 10 (e.g., Kellermann et 
al. 1989). PSS 0121+0347 is radio loud and shows the flattest 
a ox value in our sample, consistent with results from previous 
studies (e.g., Zamorani et al. 1981); the other eight objects are 
radio quiet. 



3.1. X-ray flux versus AB 1450 ( 1+Z ) magnitude 

Chandra and other exploratory X-ray observations have de- 
fined the typical fluxes and luminosities of z > 4 quasars (e.g., 
Kaspi et al. 2000; V01; Brandt et al. 2002a, 2002b; Bech- 
told et al. 2002). Figure 4 shows the observed-frame, Galactic 
absorption-corrected 0.5-2 keV flux versus AB 1450 ( 1+Z \ magni- 
tude for a compilation of z > 4 Active Galactic Nuclei (AGNs). 
The objects presented in this paper are shown as filled trian- 
gles, while open triangles and large downward-pointing arrows 
are previous Chandra observations of high-redshift quasars 
(e.g., V01; Brandt et al. 2002a), mostly taken from the 
SDSS. Circled triangles are RLQs, according to the definition 
in §3. In the following analyses and plots we exclude all but 
one (SDSS 120441.73-002149.6) of Bechtold Cycle 2 z > 4 
quasars (Bechtold et al. 2002) since they were measured in 
a different way than all other objects presented here. The 
BOSAr-detected RQQs and RLQ are plotted as filled circles 
and a circled star, respectively (Kaspi et al. 2000; V01), while 
small downward-pointing arrows indicate ROSAT upper limits, 
most of which are presented here for the first time (see Ap- 
pendix A for further information). Blazars from the Kaspi et 
al. (2000) sample are plotted as open stars. One object (the 
z = 4.01 blazar GB 1713+2148; Hook & McMahon 1998) has 
its X-ray data published here for the first time and is a w 3o 
ROSAT HRI detection. As a comparison sample, in Fig. 4 
we also plot as open circles the z = 4.45 X-ray selected RQQ 
RX J1052+5719 (Schneider et al. 1998), the z = 5.74 Broad 
Absorption Line (BAL) quasar SDSS 1044-0125 (Fan et al. 



7 h keV has been computed from the observed-frame 0.5-2 keV flux (assuming an X-ray power-law photon index of T = 2); this is required because ACIS has limited 
spectral response below 0.5 keV. Thus, the derived a OK values are actually based on the relative amount of X-ray flux in the 0.5 ( 1 + z) keV to 2( 1 + z) keV rest-frame 
band. 



Table 3 



Properties of z > 4 Quasars Observed by Chandra 



Object 




-^^1450(1+^) 


/2500 b 


log(vLv) 2 500 


M B 


Count rate 




/2keV C 


log(vLv) 2 keV 


log(L 2 -10 kev) f 


a,„ g 


R h 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


PSS 0121+0347 


3.27 


18.5 


2.22 


46.9 


-28.4 


10.30+]™ 
5 


34.7+j* 


26.4 


45.6 


45.8 


-1.51+0.06 


300.2 1 


PSS 0133+0400 


3.06 


18.2 


2.93 


47.1 


-28.7 


17.45:2 


13.4 


45.3 


45.5 


-1.67+0.06 


<4A' 


PSS 0209+0517 


4.73 


17.8 


4.24 


47.2 


-29.1 


3 64+ ' 58 


12.8±I;| 


9.8 


45.2 


45.4 


1 7 n+0.07 


< 3.0 1 


PSS 0926+3055 


1.88 


16.7 


11.70 


47.7 


-30.2 


8.67+15 


27.7^1 


21.4 


45.5 


45.7 


-1.82+0.06 


<0.3i 


PSS 0955+5940 


1.28 


18.4 


2.44 


47.0 


-28.6 


1 CQ+0.72 
1 - :>S -0.51 


5.0+f 2 


3.9 


44.8 


45.0 


-1.84+0.08 


<5.4 j 


PSS 0957+3308 


1.52 


18.2 


2.93 


47.1 


-28.7 


O Q8+0-88 
Z - VO -0.69 


9Atg 


7.3 


45.1 


45.3 


-1.77+0.07 


< LP 


PSS 1326+0743 


2.01 


17.2 


7.36 


47.4 


-29.7 


7.621UI 


24.41" 


18.6 


45.4 


45.7 


-1.77+0.06 


<0.5J 


PSS 1347+4956 


1.23 


17.4 


6.12 


47.4 


-29.7 


A 5Q+1-07 


143%* 


11.8 


45.3 


45.5 


1 O1+0.07 
_1 -° l -0.06 


0.1 k 


PSS 1443+5856 


1.50 


17.7 


4.64 


47.3 


-29.3 


I o 1+0.83 
'■ zl -0.45 


3*1™ 


3.0 


44.7 


44.9 


, qq+0.10 


<2.8 ] 





Note. — Luminosities are computed using H„=1Q km s" 1 Mpc" 1 , £2 M =0.3, and £2 A =0.7. 

"From Dickey & Lockman (1990), in units of 10 20 cirr 2 . 

b 2500 A rest-frame flux density, in units of 10~ 27 erg cm~ 2 s _1 Hz -1 . 

^Observed count rate computed in the 0.5-2 keV band, in units of 10~ 3 counts s _1 . 

d Galactic absorption-corrected flux in the observed 0.5-2 keV band, in units of 10~ 15 erg cm -2 s _1 . 

°Rest-frame 2 keV flux density, in units of 10~ 32 erg cm -2 s _1 Hz -1 . 

f 2-10 keV rest-frame luminosity corrected for Galactic absorption, in units of erg s _1 . 

s Errors have been computed following the "numerical method" described in § 1.7.3 of Lyons (1991); both the statistical uncertainties on the X-ray count rates and the effects of the 
observed ranges of the X-ray and optical continuum shapes have been taken into account (see text for details; see also §3 in V01). 

h Radio-loudness parameter, defined as R = f 5 ghz// 4400 ^ (rest frame) (e.g., Kellermann et al. 1989). The rest-frame 5 GHz flux density is computed from the observed 1 .4 GHz 
flux density assuming a radio power-law slope of a = -0.8, with f v <x v a . 

n.4 GHz flux density from NVSS (Condon et al. 1998). 

J 1.4 GHz flux density from FIRST (Becker et al. 1995). 

k 1.4 GHz Very Large Array flux density from Carilli et al. (2001a). 
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2000; Brandt et al. 2001b), and the z = 4.42 Seyfert galaxy 
VLA J1236+6213 (Waddington et al. 1999; Brandt et al. 
2001c). The three filled squares show two quasars and a Seyfert 
galaxy serendipitously discovered by Chandra (Silverman et al. 
2002; Barger et al. 2002) at z = 4.93, z = 5.18, and z = 4.13. 
The slanted lines show a ox spectral indices of — 1 .5 and — 1 .8 at 
z = 4.2, which is the average redshift for the sample discussed 
in this paper. 




SDSS J1044-0125 



CDF-N z=4.13- CDF-N z = 5.^' 



VLA J1236+62KU " 

io ": pr a 

» y ' 

10-" U 1 1 1 — 

26 24 



I 

22 



20 

magnitude 





18 



FIG. 4. — Observed-frame, Galactic absorption-corrected 0.5-2 keV flux 
versus ABinj ( i+z ^ magnitude for z > 4 AGNs. The objects presented in 
this paper are plotted as filled triangles, while the open triangles and large 
downward-pointing arrows indicate previous Chandra observations of high- 
redshift quasars (V01; Brandt et al. 2002a). Circled triangles are RLQs. The 
ROSAT-detected RQQs and RLQ are plotted as filled circles and a circled star, 
respectively (Kaspi et al. 2000; V01), while small downward-pointing arrows 
indicate ROSAT upper limits (most of the upper limits are presented here for 
the first time; see Appendix A for details). Blazars from the Kaspi et al. (2000) 
sample are plotted as open stars (one object is published here for the first time; 
see Appendix A). For comparison, we have plotted as open circles the z = 4.45 
X-ray selected RQQ RX J1052+5719 (Schneider et al. 1998), the z = 5.74 
RQQ SDSS J1044-0125 (Fan et al. 2000; Brandt et al. 2001b), and the 
Z = 4.42 Seyfert galaxy VLA J1236+6213 (Waddington et al. 1999; Brandt et 
al. 2001c). The three filled squares indicate z > 4 AGNs serendipitously dis- 
covered by Chandra (Silverman et al. 2002; Barger et al. 2002). The slanted 
lines show a ox =-1.5 and a ox =-1.8 loci at z = 4.2 (the average redshift of the 
present sample), assuming the same X-ray and optical spectral shapes used in 
the text. 

As shown in Fig. 4, quasars at z > 4 are generally faint X-ray 
sources, with 0.5-2 keV fluxes <4x 10~ 14 ergcm _2 s _1 . 
The PSS quasars presented in this paper are among the X-ray 
brightest known and allow study of a region of the luminosity- 
redshift parameter space different from that of most of the 
SDSS quasars observed in X-rays to date (e.g., V01; Brandt 
et al. 2002a), which are usually optically less luminous and 
located at higher redshift. Given the X-ray fluxes found from 
Chandra snapshot observations, it is clear that PSS and BRI 8 
quasars represent the best known RQQs at z > 4 for longer 
XMM-Newton spectroscopic observations. Individual XMM- 
Newton spectra, coupled with stacked Chandra spectra, can 
provide tighter spectral constraints and possibly allow searches 



for spectral features such as iron Ka lines. 

In Fig. 4 there is a correlation between the soft X-ray flux 
(corrected for Galactic absorption) and AB 14 5 ( 1+Z ) magni- 
tude. Given the presence of upper limits, to study this cor- 
relation we have used the ASURV software package Rev 1.2 
(La Valley, Isobe, & Feigelson 1992). To evaluate the signif- 
icance of this correlation, we used several methods available 
in ASURV, namely the Cox regression proportional hazard 
(Cox 1972), the generalized Kendall's x (Brown, Hollan- 
der, & Korwar 1974), and the Spearman's p models. Us- 
ing the Chandra RQQs, the correlation is significant at the 
99.5-99.8% and 99.99% levels when all the PSS/BRI and 
PSS/BRI/SDSS quasars are taken into account, respectively. 
In the latter case, the best-fit relationship is parameterized by 
fx = [-(6.9 ±1.0) AB 1450(1+Z) + (137 ± 19)] x IO" 15 erg cm- 
Note that the scatter around the best-fit relationship is signifi- 
cant (a factor of w 3). The BALQSOs have been excluded in 
these correlation analyses because they are known to be signif- 
icantly absorbed in the X-ray band (e.g., Brandt, Laor, & Wills 
2000; Gallagher et al. 2001, 2002; Green et al. 2001; see also 
V01 for discussion about the BALQSOs previously observed 
and undetected by Chandra at z > 4). Including also the RQQs 
with ROSAT detections or tight upper limits (i.e., a ox < —1.5) 
provides a significance > 99.99% for the above correlation. 
The observed correlation indicates that the optical and X-ray 
emission are likely to be powered by the same engine, namely 
accretion onto a supermassive black hole. The large 2-10 keV 
rest-frame luminosities of the PSS objects investigated in this 
paper (« io 449 ~ 45 - 8 erg s _1 ) suggest that the X-ray emission 
is mostly nuclear and not related to the starburst component 
whose existence has been suggested by recent observations at 
millimeter and sub-millimeter wavelengths (e.g., Omont et al. 
2001; Isaak et al. 2002; see also Elvis, Marengo, & Karovska 
2002 for a different interpretation). 



3.2. a ox results 

The average a ox for the eight RQQs of the present PSS 
sample is (a ox )=— 1.81 ±0.03 (the quoted errors represent 
the standard deviation of the mean), with the BALQSO 
PSS 1443+5856 (see §4) being the quasar with the steepest 
a ox in our sample (a ox =— 1.99^g'gg). Even after removing the 
BALQSO, the average a ox is -i.78±0.02, still considerably 
steeper than that obtained for samples of lower-luminosity, lo- 
cal quasars [e.g., the Bright Quasar Survey (BQS; Schmidt 
& Green 1983) quasars at z < 0.5 have (a ox }=-1.56±0.02; 
see Brandt et al. 2000 and V01]. The same applies 
when all the PSS/BRI ((a ox )=-1.78±0.03) or the SDSS 
((oc ox )=— 1.75±0.03) z > 4 quasars observed by Chandra are 
taken into account. These average a ox values are slightly flat- 
ter (typically by 0.02-0.03) when the BALQSOs are excluded 
from the analysis. One possible explanation for the more nega- 
tive a ox values seen at z > 4 is the presence of X-ray absorption; 
however, we find no evidence of X-ray absorption either from 
the individual quasars' spectra (see Table 2) or from the joint 
spectral fitting discussed in §3.3. Another possibility is that a ox 
depends upon 2500 A rest-frame luminosity (e.g., Avni, Wor- 
rall, & Morgan 1995) or redshift (e.g., Bechtold et al. 2002). 
We have searched for significant correlations of a ox with either 



8 BRI quasars have been selected from the B/R/I survey (Irwin, McMahon, & Hazard 1991) and then confirmed spectroscopically (e.g., Storrie-Lombardi et al. 1996, 
2001). Their selection criteria and optical magnitude range are similar to those of PSS quasars. See V01 and Bechtold et al. (2002) for the X-ray properties of BRI 
quasars in Chandra observations. 



X-rays from the most luminous optically selected z > 4 quasars 
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FIG. 5. — (a) a ox versus log ^ 25f|0 ^ and (b) redshift for z > 4 optically 
selected RQQs. Symbols are the same as in Fig. 4. 



Z > 4 using the statistical methods available in ASURV (see 
§3.1) but found none. Given the small spread in UV luminosi- 
ties and redshifts of the objects used in this analysis, we cannot 
address the question of a ox changes with rest-frame 2500 A lu- 
minosity or redshift effectively. However, the analysis of SDSS 
quasars in the redshift range z « 0.2-6.3 using partial correla- 
tion analysis (Vignali, Brandt, & Schneider 2002) shows that 
a ox is mainly dependent upon 2500 A luminosity. 

3.3. Joint spectral fit 

To provide tighter constraints on the average X-ray spectral 
properties of z > 4 quasars, we have performed a joint spectral 
analysis of our target quasars. Source counts have been ex- 
tracted from 2" radius circular apertures centered on the X-ray 
position of each quasar. The background has been taken from 
annuli centered on the sources, avoiding the presence of nearby 
faint X-ray sources. To account for the recently discovered 
quantum efficiency decay of AGS 9 at low energies, possibly 
caused by molecular contamination of the ACIS filters, we have 
applied a time-dependent correction to the ACIS quantum effi- 
ciency generating a new ancillary response file (ARF) for each 
source. 10 The quantum efficiency degradation is not severe at 
energies above w 0.7 keV (« 10% at 1 keV), where most of 
the detected signal lies, but, if left uncorrected, it might affect 
the results at lower energies. In particular, the reduction in ef- 
ficiency could produce unreliable column densities if not taken 
into account properly. 

Only two objects are characterized by low counting statis- 

9 See http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/index.html. 

10 See http://www.astro.psu.edu/users/chartas/xcontdir/xcont.html. 



tics (« 10 counts); all the others have w 20-70 counts in the 
0.5-8 keV band. The sample does not seem to be biased by 
the presence of either a few high signal-to-noise ratio objects 
(see Table 2) or the RLQ PSS 0121+0347, whose X-ray pho- 
ton index (derived from the band ratio analysis; see Table 2) 
is consistent with those of the other quasars. Spectral analy- 
sis was carried out with XSPEC. The joint spectral fitting of 
the nine individual unbinned quasar spectra was carried out 
using the Cash statistic (Cash 1979). The Cash statistic is 
particularly well suited to low-count sources (e.g., Nousek & 
Shue 1989), and, since it is applied to the unbinned data, it 
allows us to retain all the spectral information. Furthermore, 
we can apply to each source its own Galactic absorption and 
redshift. The resulting joint spectrum, although derived from 
only « 340 counts, is reasonably good due to the extremely 
low background of Chandra (there are only sa 2. 1 background 
counts in all source cells), and it allows spectral analysis in the 
w 2-30 keV rest-frame band. Joint fitting with a power-law 
model (leaving the normalizations free to vary) and Galactic 
absorption provides a good fit (as judged by the small data- 
to-model residuals and checked with 10000 Monte-Carlo sim- 
ulations) with a photon index of r=1.98±0.16 (errors are at 
the 90% confidence level; A% 2 =2.71, Avni 1976). This pho- 
ton index is consistent with those obtained for lower-redshift 
samples of RQQs in the 2-10 keV band (e.g., George et al. 
2000; Mineo et al. 2000; Reeves & Turner 2000) and with 
that of the BALQSO APM 08279+5255 at z = 3.91 (Chartas 
et al. 2002; F=U2^q'^). The photon index is also consis- 
tent with the average F determined from the band ratios (see 
Table 2) using ASURV ((r)=2.07±0.13). The combined spec- 
trum (used here only for presentation purposes), fitted with the 
power-law model described above, is shown in Fig. 6. Bech- 
told et al. (2002) reported an X-ray spectral flattening with red- 
shift for quasars, with photon indices in the range F « 1.2-1.6 
((r}=1.48±0.13) at z « 3.7-6.3. The present data, obtained for 
an optically brighter sample of quasars, are not consistent with 
such claims. 




0.5 1 1.5 
N (10 2s cuT 3 ) 



0.5 1 2 5 

Energy (keV) 

FIG. 6. — Combined quasar spectrum (used only for presentation purposes) 
fitted with the best-fit power-law model and Galactic absorption (see §3.3 for 
details). The 68, 90, and 99% confidence regions for the photon index and in- 
trinsic column density obtained using the Cash statistic are shown in the insert. 
Data-to-model residuals are shown in the bottom panel (in units of a). 
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Furthermore, we have constrained the presence of neutral in- 
trinsic absorption. In the absorption column we have assumed 
solar abundances, although there are indications of supersolar 
abundances of heavy elements in high-redshift quasar nuclei 
(e.g., Hamann & Ferland 1999; Constantin et al. 2002). Dou- 
bling the abundances in the fit gives a reduction in the column 
density of a factor of w 2. There is no evidence for absorption 
in excess of the Galactic column density, as shown in the in- 
sert in Fig. 6, where the 68, 90, and 99% confidence regions 
for the photon index and intrinsic column density are plotted. 
Although it is not possible to rule out substantial X-ray absorp- 
tion in the individual spectrum of any one quasar, on average 
it appears that any neutral absorption has a column density of 
Mi ~ 8.8 x 10 21 crrr 2 (at 90% confidence). Note that this av- 
erage constraint is more generally applicable than if we had just 
a « 340-count spectrum of a single z > 4 quasar, since any one 
quasar might always be exceptional. If our targets had intrin- 
sic absorption of Na s=s (2-5) x 10 22 cm~ 2 , as recently found in 
two RQQs atz«2 using ASCA data (Reeves & Turner 2000), 
we would be able to detect it given the upper limit on the col- 
umn density we derived. We have also investigated the possible 
presence of ionized absorption using the model ABSORI in 
XSPEC. While, of course, this is more difficult to constrain in 
general, our data do not provide any suggestion of ionized ab- 
sorption. Furthermore, considering narrow iron Ka lines in the 
rest-frame energy range 6.4-6.97 keV, the typical upper limit 
on the equivalent width (EW) is « 420 eV. 

3.4. Companion objects 

Schwartz (2002b) tentatively associated an X-ray source 23" 
from the z = 5.99 quasar SDSS 1306+0356 with an X-ray jet, 
although neither the source nor SDSS 1306+0356 have any ra- 
dio emission in the FIRST catalog (Becker et al. 1995). To 
date, no quasar with X-ray emission undoubtedly associated 
with jets is radio quiet, all being powerful radio sources (e.g., 
Harris & Krawczynski 2002). Moreover, the jets are char- 
acterized by substantial radio emission (e.g., Siemiginowska 
et al. 2002; Brunetti et al. 2002 and references therein) 
even after correction for boosting. If the X-ray emission from 
jets is produced via the Compton scattering of cosmic mi- 
crowave background (CMB) photons by relativistic electrons, 
then the enhancement of the CMB energy density with red- 
shift compensates for the decrease of surface brightness (see 
Schwartz 2002a for details; also see Rees & Setti 1968). Re- 
cently, Ivanov (2002) detected the /-band counterpart of the 
X-ray source 23" from SDSS 1306+0356. The elongation of 
the /-band counterpart suggests the presence of a foreground 
galaxy, rather than UV/optical emission from a jet associated 
with SDSS 1306+0356, but no redshift for the counterpart has 
yet been measured. 

In the present sample of nine objects, only one quasar, 
PSS 0121+0347, is radio loud (see §3). Nevertheless, we have 
searched for possible companions/jets for all the objects in our 
sample over a region of w 100"x 100" centered on the quasars 
(the same shown in the images of Fig. 2). At z ~ 4^1.5, this 
corresponds to a linear scale of w 320 kpc at each side of the 
quasar. We found some faint objects (about one per field); none 
of these is X-ray bright and close enough to our targets to con- 
taminate their X-ray emission. Most of the serendipitous X-ray 
sources lack counterparts in the DPOSS images, and none of 



them has associated radio emission down to a 1 .4 GHz flux den- 
sity of « 0.5-1.5 mJy (3a). 

As an observational check, we compared the 0.5-2 keV cu- 
mulative number counts of sources detected in the current 
high-redshift quasar fields with the number counts obtained 
from the Chandra Deep Field-North survey (Brandt et al. 
2001a) and the ROSAT deep survey of the Lockman Hole 
(Hasinger et al. 1998); the high-redshift quasars themselves, of 
course, were excluded. The surface density of X-ray sources 
surrounding the PSS quasars at the flux limit we reached 
(« 10~ 15 erg cm" 2 s" 1 ) is N(> S)=1224^ 2 deg" 2 , consis- 
tent with the densities obtained in deep-field observations. Thus 
we find no evidence for an excess of companions/jets associated 
with our quasars. 

3.5. X-ray extension of the quasars 

To constrain the presence of either gravitational lensing (e.g., 
Barkana & Loeb 2000; Wyithe & Loeb 2002; Comerford, 
Haiman, & Schaye 2002) or X-ray jets close to the sources (e.g., 
Schwartz 2002b), we performed an analysis of X-ray extension 
for all the quasars in the present sample. At z « 4-5, the angu- 
lar separation between the images of a multiply imaged source 
should often be w 1-2", depending on the adopted dark mat- 
ter model (Barkana & Loeb 2000), and Chandra has been able 
to detect and separate lensed images up to z « 3.9 (Chartas et 
al. 2002). For each quasar, we created a point spread func- 
tion (PSF) at w 1.5 keV at the source position (using the tool 
MKPSF) normalized by the effective number of source counts. 
Then we compared the radial profile of this PSF with that of the 
source, finding good agreement in the radial counts distribution. 
For one quasar (PSS 0926+3055) we performed a further check 
given the presence in the soft band of three counts "pointing" 
in the North- West direction within a radius of 3" of the source, 
although there is no indication of extension in the optical im- 
age. We deconvolved the X-ray source image with the PSF and 
found that the counts distribution still seems to be in agreement 
with the source being point-like. However, further investigation 
with a longer Chandra observation is needed to completely rule 
out the possibility of PSS 0926+3055 being slightly extended. 

3.6. X-ray variability of the quasars 

Although the present observations are short (the observation 
length in the rest frame is w 15 minutes), we searched for 
X-ray variability of our targets. There are tentative reports that 
quasar X-ray variability increases with redshift (e.g., Manners, 
Almaini, & Lawrence 2002). We analyzed the photon arrival 
times in the 0.5-8 keV band using the Kolmogorov-Smirnov 
(KS) test. No significant variability was detected (the KS prob- 
abilities that the sources are variable range from « 6% to 83%). 
Our constraints on the amplitude of variability depend upon 
the number of counts detected as well as the putative variabil- 
ity timescale. Systematic variability sustained over a period of 
w 5-6 ks (in the observed frame) must be less than w 30-90%. 

4. OPTICAL ANALYSIS 

Following the method described in §2.2, two quasars, 
PSS 0926+3055 and PSS 0955+5940, have been found to vary 
in the optical band on a timescale of w 6 years (in the observed 



Table 4 
Emission-Line Measurements 





Lya(1216 A) 


NV(124()A 


I 




01(1304 A) 




Si IV-0 IV] (1394-1402 A) 




CIV (1549 A) 




Object Name 


X' EW b 


X* EW b 


z 


X' 


EW b FWHM C 


z 


X" EW b FWHM C z 


X a 


EW b FWHM C 


(z) 



PSS 0121+0347 
PSS 0133+0400 
PSS 0209+0517 
PSS 0926+3055 
PSS 0955+5940 
PSS 0957+3308 
PSS 1326+0743 
PSS 1347+4956 
PSS 1443+5856 



6232.3 
6266.1 
6326.1 
6309.0 
6503.5 
6386.9 
6223.6 
6760.6 
6465.0 



81.1 
25.9 
33.5 
44.0 
285.0 
22.1 
26.2 
8.7 
69.5 



6391.7 2.3 4.152 



6306.3 4.2 4.083 



6674.4 


4.8 


1020 


4.112 


6734.8 


10.3 


6270 


4.158 


6802.6 


1.3 


1460 


4.210 


6979.2 


3.9 


2270 


4.345 


6829.5 


1.2 


2150 


4.231 


6670.6 


1.6 


2260 


4.109 



7218.4 


5.3 


2550 


4.150 


7216.8 


1.1 


3350 


4.149 


7261.6 


8.2 


6620 


4.181 


7471.1 


4.0 


2730 


4.330 


7267.2 


9.4 


6340 


4.185 


7122.7 


10.0 


7130 


4.082 


7708.4 


4.1 


5630 


4.500 


7369.2 


9.9 


6760 


4.258 



7963.6 
7999.7 
7945.1 
8028.1 
8270.0 
8020.7 
7905.4 
8535.3 
8161.8 



33.5 
32.6 
12.0 
18.5 
21.6 
12.1 
21.2 
15.9 
19.0 



3810 
7370 

12940 
5880 
4220 
9500 

10920 
7670 
7980 



4.142 
4.165 
4.130 
4.184 
4.340 
4.179 
4.105 
4.511 
4.270 



4.126 
4.155 
4.143 
4.191 
4.339 
4.198 
4.089 
4.506 
4.264 



Note. — The redshift for each quasar was determined by averaging the emission lines (excluding Lya for all objects). The last column reports the weighted average for the redshift. The FWHMs of Lya and N V are not 
reported because of the uncertainties in their measurements (see, e.g., Schneider et al. 1991). 

"Observed-frame wavelength (A). 

b Rest-frame equivalent width (A). 

TWHM of the line (km s~'), corrected for the instrumental resolution. 
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frame). The HET /?-band magnitude for the former is 0.3 mag 
brighter than the value reported in the Djorgovski list (taking 
into account the different filter transmission and bandpass; see 
Footnote 4), while the latter appears 0.5 mag fainter. These are 
the only objects for which significant long-term optical vari- 
ability has been found. 

The principal emission-line measurements are presented in 
Table 4. The centroids of the emission lines were determined 
with Gaussian fits. Redshifts have been derived from HET 
spectra averaging the emission-line redshifts (excluding Lya, 
whose measurements are affected by many uncertainties; see, 
e.g., Schneider, Schmidt, & Gunn 1991). Typical uncertain- 
ties are w 0.01 for non-BALQSO redshift measurements, and 
w 10-15% and w 5-10% for FWHM and EW measurements, 
respectively. We have compared the C IV emission-line equiva- 
lent widths of our objects with those derived for the Constantin 
et al. (2002) sample of optically luminous z > 4 PSS and 
BRI quasars; we found good agreement. One of our quasars, 
PSS 1443+5856, shows broad C IV and Si IV absorption fea- 
tures in its optical spectrum (see Fig. 3) with a blueshift of 
w 10000 km s -1 with respect to the systemic redshift of the 
quasar. The C IV absorption-line EW is « 14 A (rest frame), 
which is close to the value expected from the a 0X -C IV EW 
anti-correlation found for low-redshift quasars (see Fig. 4 of 
Brandt et al. 2000). The likely existence of this anti-correlation 
at z > 4 has already been verified (e.g., Maiolino et al. 2001; 
Goodrich et al. 2001; also see V01). Given the correlation 
between AB l4 _ 50 ^ 1+z - ) and soft X-ray flux shown in Fig. 4, the 
X-ray detection of PSS 1443+5856 is not surprising. Since 
this quasar is 1.5-1.7 magnitudes brighter than the previously 
X-ray undetected z > 4 BALQSOs (V01), a higher X-ray flux 
is expected. 

5. SUMMARY AND IMPLICATIONS 

We have presented the Chandra X-ray detections of nine 
high-redshift (zw 4.1-4.5) quasars. These quasars were tar- 
geted because they are among the most optically luminous 
objects in the Universe and are likely to comprise a signif- 
icant fraction of the most luminous z > 4 optically selected 
objects. Their SEDs are characterized by steeper (more neg- 
ative) optical-to-X-ray spectral indices ((a ox }=— 1.81+0.03) 
compared to samples of optically selected quasars at low or in- 
termediate redshift. Given their large mean 2500 A luminosity 
density (« 1.6x 10 32 erg s -1 Hz -1 ), their steep a ox values may 
be explained in the context of the a ox -ultraviolet luminosity 
anti-correlation found in other work (e.g., Avni et al. 1995; Vi- 
gnali et al. 2002). We also confirm the presence of a significant 
correlation between the optical magnitude and soft X-ray flux 
previously found by V01 for a more heterogeneous sample of 
z > 4 quasars. Joint spectral fitting in the 2-30 keV rest- 
frame band indicates that optically selected, luminous z > 4 
quasars are characterized, on average, by a power-law model 
with a photon index of w 2.0+0.2, similar to the majority of 
lower-redshift quasars. There is no evidence for intrinsic ab- 
sorption (Na ;$ 8.8 x 10 21 crrr 2 at the 90% confidence level). 

We also have presented near-simultaneous optical spectra ob- 
tained with the HET. Optical spectra for seven of these quasars 
are published here for the first time. One of the quasars, 
PSS 1443+5856, shows broad absorption features in C IV and 
Si IV. 

From a more general perspective, our results show that 
z ~ 4.1-4.5 quasars and local quasars have reasonably simi- 



lar X-ray and broad-band spectra (once luminosity effects are 
taken into account). The small-scale X-ray emission regions of 
quasars do not appear to be sensitive to the substantial large- 
scale environmental differences between z ~ 4 and z~0. The 
lack of any strong X-ray spectral changes is particularly rele- 
vant since X-ray emission originates from the compact region 
where accretion, and thus black hole growth, occurs. Our data 
do not provide any hints for different accretion/growth mech- 
anisms (via, e.g., accretion-disk instabilities or "trapping ra- 
dius" effects) between the most luminous high-redshift quasars 
and local quasars. Furthermore, our data suggest that (1) min- 
imum black hole masses for high-redshift quasars can be es- 
timated reasonably well via the Eddington limit by adopting 
the bolometric correction used for comparably luminous low- 
redshift quasars, and (2) calculations of the effects of early 
quasars upon the intergalactic medium can, at least to first or- 
der, adopt the spectral energy distributions of comparably lumi- 
nous low-redshift quasars. Finally, our data support the dogma 
that X-ray emission is a universal property of quasars at all red- 
shifts. Quasars at even higher redshifts (z m 5-15) will likely 
be luminous X-ray sources that can be detected in future X-ray 
surveys. 



X-rays from the most luminous optically selected z > 4 quasars 
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FIG. 7. — (a) Simulated 40 ks spectrum and (b) 68, 90, and 99% confidence 
regions for the He-like iron Ka emission-line energy and intensity using the 
calorimeter and the grating onboard Constellation-X. We have used the ob- 
served Chandra 0.5-2 keV flux of ss 2.8 x 10~ 14 erg cm~ 2 s _1 and assumed 
a r = 2.0 power-law spectrum with Galactic absorption. We also included a 
narrow iron Ka emission line with a rest-frame equivalent width of 70 eV (see 
§5 for details), (c) 68, 90, and 99% confidence regions for intrinsic column 
density and photon index. The "X" marks the input parameters of the simula- 
tion, while the "+" indicates the best-fit result. 



6. FUTURE WORK 

We have recently been awarded 1 1 additional Chandra obser- 
vations of z > 4 optically luminous quasars. These observations 
will provide tighter constraints on the X-ray properties of high- 
redshift quasars, thus allowing for more feasible planning of 
X-ray observations with next-generation, large-area X-ray tele- 
scopes. In this regard, we have simulated a 40 ks Constellation- 
X il spectrum (with the calorimeter and the first-order grat- 
ing; see Fig. 7a) of the X-ray brightest RQQ of our sample, 
PSS 0926+3055 at z = 4.19, assuming a power-law spectrum 
with r=2.0 and Galactic absorption. A narrow (o=10 eV), He- 
like (at 6.7 keV) iron Ka emission line with a rest-frame equiv- 
alent width of « 70 eV has also been included in the model. At 
high X-ray luminosities (L2-10 kev>10 44 ~ 45 erg s _1 ) the iron 
Ka line is expected to be ionized and fairly weak, with most of 
the line flux coming from its core (the "X-ray Baldwin effect" 
discussed in Nandra et al. 1997; see also Iwasawa & Taniguchi 
1993). The 68, 90, and 99% confidence regions for the line 
energy and intensity are plotted in Fig. 7b. Constellation-X 
would gather w 11300 counts (fa 8500 from the calorime- 
ter and ;=y 2800 from the first-order grating) in the observed 
0.25-10 keV band. These would ensure a good accuracy in the 
measurements of the X-ray parameters; in particular, the pho- 
ton index, line energy, and line EW would be measured with 
accuracies of w 1.2%, 0.1%, and 29%, respectively. We are 
able to detect ionized iron Ka lines down to rest-frame EWs of 
~50 eV. Due to its low-energy coverage, Constellation-X can 
also place significant constraints on the intrinsic column den- 
sity down to Nh ~ 2 x 10 21 cm~ 2 (see Fig. 7c). 
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APPENDIX 

NEW ROSAT UPPER LIMITS FOR Z > 4 QUASARS AND A TENTATIVE DETECTION OF GB 1713+2148 

The Djorgovski list of z > 4 quasars (see Footnote 4) has been cross correlated with ROSAT PSPC and HRI archival observations. 
In the 0.5-2 keV energy band, we found 45 new upper limits [34 RQQs, five RLQs, four quasars with loose radio loudness upper 
limits (R £ 15), and two quasars without FIRST/NVSS coverage] and a tentative w 3o detection for the blazar GB 1713+2148 (in 
an HRI observation). The 3o upper limits (reported in Table Al) have been derived in a similar way to those shown in Table 4 of 
V01 using the SOSTA command in the XIMAGE package (Giommi et al. 1992). For the 3a detection of GB 1713+2148, we cross 
checked the results obtained with the maximum likelihood method available with the MIDAS/EXS AS package (Zimmermann et al. 
1998) with the XIMAGE results. GB 1713+2148 is unusually X-ray faint compared to the other blazars, despite its large radio flux 
density (« 450 mJy in the NVSS; « 24400). 

For SDSS quasars, the AB 145 o(i +z ) magnitudes have been derived from the /-band filter using either the empirical relationship 
ABi45o(i +z )=z' — 0.2 or the method described in Vignali et al. (2002) for the SDSS quasars in the Early Data Release (Schneider et al. 
2002). For the other quasars, the AB l450 ^ l+z - > magnitudes have been derived either from the /?-band magnitude using the empirical 
relationship AB 145Q ^ 1+Z ^=R— 0.684 z + 3.10 or from the r-band magnitude using AB 1450 ( 1+Z )=r— 0.684 z + 2.75. Unless otherwise 
stated, the R- and r-band magnitudes have been obtained from the Djorgovski list (see Footnote 4). All of the above relationships 
provide reliable AB 14 5 ( 1+Z ) estimates (within w 0. 1-0.2 magnitudes) in the redshift range under consideration. 
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Table A 1 

Properties of z > 4 Quasars in ROSAT Fields 



Object 


z 


AB 1450(l+z) 


M B 




ct 


£ 




ROSAT RQQs 


PSS 0014+3032 


4.47 


18.5 


-28.6 


<4.44 


< -1 


.45 


BRI 0046-2458 


4.15 


18.9 


-28.1 


< 3.97 


< -1 


.43 


SDSS 0127-0045 


4.06 


17.9 


-29.0 


< 3.53 


< -1 


.60 


PSS 0134+3307 


4.53 


18.5 


-28.6 


< 13.3 


< -1 


.27 


PSS 0223+2325 


4.27 


18.8 


-28.2 


< 12.3 


< -1 


.25 


SDSS 0239-0810 


4.00 


19.4 


-27.4 


< 1.90 


< -1 


.47 


SDSS 0252+0031 


4.10 


19.7 


-27.2 


< 11.8 


< -1 


.13 


BR 0300-0207 


4.25 


18.4 


-28.6 


< 8.46 


< -1 


.37 


BR 0353-3820 


4.55 


17.6 


-29.5 


< 11.7 


< -1 


.43 


FIRST 0747+2739 b 


4.11 


18.2 


-28.7 


< 11.3 


< -1 


.36 


PSS 0808+5215 


4.44 


18.5 


-28.5 


<4.58 


< -1 


.45 


SDSS 0810+4603 


4.07 


18.2 


-28.7 


< 7.03 


< -1 


.55 


SDSS 0832+5303 


4.02 


19.6 


-27.3 


<22.0 


< -1 


.08 


PSS 0852+5045 


4.20 


18.0 


-29.0 


<7.24 


< -1 


.47 


SDSS 0904+5350 


4.25 


19.1 


-27.9 


< 3.90 


< -1 


.40 


PSS 1026+3828 


4.18 


18.8 


-28.1 


< 3.20 


< -1 


.47 


SDSS 1034-0027 


4.38 


19.8 


-27.2 


<7.48 


< -1 


.17 


SDSS 1040-0015 


4.32 


18.7 


-28.3 


<8.28 


< -1 


.32 


SDSS 1048-0028 


4.00 


19.2 


-27.7 


<4.66 


< -1 


.35 


SDSS 1052-0006 


4.13 


19.6 


-27.3 


<3.33 


< -1 


.34 


SDSS 1056+0032 


4.02 


19.8 


-27.1 


<9.52 


< -1 


.15 


SDSS 1059+0104 


4.06 


19.4 


-27.6 


< 8.61 


< -1 


.23 


PSS 1140+6205 


4.51 


18.4 


-28.7 


<2.86 


< -1 


.54 


PSS 1226+0950 


4.34 


18.6 


-28.4 


<7.04 


< -1 


.37 


PSS 1248+3110 


4.32 


19.1 


-27.9 


<4.60 


< -1 


.37 


PSS 1315+2924 


4.18 


18.4 


-28.6 


< 1.66 


< -1 


.65 


PSS 1339+5154 


4.08 


19.0 


-27.9 


< 1.21 


< -1 


.61 


PSS 1401+4111 


4.01 


18.6 


-28.2 


< 13.2 


< -1 


.27 


PSS 1418+4449 


4.32 


18.5 


-28.5 


<3.69 


< -1 


.48 


PSS 1430+2828 


4.30 


19.5 


-27.5 


<5.38 


< -1 


.28 


PSS 1543+3417 


4.41 


18.5 


-28.5 


< 13.7 


< -1 


.27 


PSS 1555+2003 


4.22 


19.1 


-27.8 


<4.26 


< -1 


.38 


PSS 1745+6846 


4.13 


19.1 


-27.9 


<0.87 


< -1 


.65 


PSS 2155+1358 C 


4.26 


18.2 


-28.8 


<7.73 


< -1 


.42 



ROSAT RLQs 



SDSS 0300+0032 


4.19 


19.9 


-27.0 


< 


5.19 


< 


-1.22 


PSS 0439-0207 


4.40 


18.8 


-28.2 


< 


5.55 


< - 


-1.37 


SDSS0913+5919 d 


5.11 


20.3 


-27.0 


< 


8.70 


< 


-1.05 


CLASS 0915+068 


4.19 


18.9 


-28.0 


< 


7.08 


< 


-1.32 


CLASS 1322+116 


4.40 


18.9 


-28.1 


< 


7.47 


< 


-1.31 



ROSAT quasars with loose radio loudness upper limits 



SDSS 0338-RD657 c 


4.96 


21.1 -26.1 


< 8.63 


< 


-0.93 


SDSS 0405+0059 


4.05 


19.8 -27.1 


< 10.0 


< 


-1.13 


SDSS 1708+6022 


4.35 


19.8 -27.2 


<2.60 


< 


-1.34 


SDSS 1710+5923 


4.47 


19.7 -27.4 


<6.00 


< 


-1.23 


ROSAT quasars without FIRST/NVSS coverage 


BR0413-4405 C 


4.07 


19.2 -27.7 


<25.1 


< 


-1.07 


BR0418-5723 C 


4.46 


18.8 -28.3 


<21.6 


< 


-1.15 


ROSAT Blazar 


GB 1713+2148 


4.01 


21.4 -25.5 


2.81 




-1.11 



Note. — The division between RQQs and RLQs is based on the radio-loudness 
parameter (defined in §3). 

"Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 
10 -14 erg cm" 2 s _1 . 

b AB 145 Q( 1+z ) magnitude derived from the R-band magnitude published in 
Richards et al. (2002); also the redshift is from Richards et al. (2002). 

C AB, 450 ( 1+Z ) magnitude derived from the R-band magnitude published in Per- 
oux et al. (2001); also the redshift for BR 0418-5723 has been obtained from 
Perouxetal. (2001). 

d Also observed by Chandra (C. Vignali et al., in preparation). 

e Also observed by XMM-ffmi'toii (C. Vignali et al., in preparation). 



